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Abstract 



We study the problem of power-efficient routing for multihop wireless ad hoc sensor networks. The guid- 
ing insight of our work is that unlike an ad hoc wireless network, a wireless ad hoc sensor network does 
not require full connectivity among the nodes. As long as the sensing region is well covered by connected 
nodes, the network can perform its task. We consider two kinds of geometric random graphs as base inter- 
connection structures: unit disk graphs UDG(2, A) and fc-nearest-neighbor graphs NN(2, k) built on points 
generated by a Poisson point process of density A in I 2 . We provide subgraph constructions for these two 
models UDG-SENS(2, A) and NN-SENS(2, k) and show that there are values A s and k s above which these 
' constructions have the following good properties: (i) they are sparse; (ii) they are power-efficient in the sense 

that the graph distance is no more than a constant times the Euclidean distance between any pair of points; 
(iii) they cover the space well; (iv) the subgraphs can be set up easily using local information at each node. 
Our analyses proceed by coupling the random graph constructions in R 2 with a site percolation process in 
Z 2 and using the properties of the latter to derive properties of the former. An important consequence of our 
constructions is that they provide new upper bounds for the critical values of the parameters A and k for the 
, models UDG(2, A) and NN(2, k). We also describe a simple local algorithm requiring only location informa- 

■ tion (from a GPS for example) and communication with immediate neighbors for setting up the subnetworks 

^ 1 UDG-SENS(2, A) and NN(2, k) and for routing packets on them. 

O . 

1 Introduction 

Multihop sensor networks, where nodes act not only to sense but also to relay information, have proven advantages 
in terms of energy efficiency over single hop sensor networks Not only is sensor-to-sensor communication 
5^ ■ useful for necessary tasks like time synchronization 6] , for certain kinds of collaborative sensing functions like 
target tracking [23| sensor-to-sensor communication is essential. The question of how to achieve connectivity 
arises here just as it does in ad hoc wireless networks with one crucial difference: 

It is not necessary that every sensor be part of a connected network. It is only necessary that the 
density of connected sensors is high enough to perform the sensing function. 

In other words, even if some sensors are wasted in the sense that the data they sense cannot be relayed, it does 
not matter as long as the area being sensed is well covered with useful sensors which are part of multihop network 
that can relay data. The difference from other ad hoc wireless networks is that each node expects connectivity 
as a service provided to it, while in the WASN individual nodes are not important, the overall task is. In this 
paper we follow this critical insight to propose sparse easy-to-compute power-efficient constructions for multihop 
WASNs. 

We consider two different types of geometric random graphs as the base interconnection structures. The 
nodes are modelled by a Poisson point process with density A on the plane, R 2 . The interconnections between 
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these points are modeled in two ways: 1) Unit Disk Graphs in which there is an edge between two nodes if the 
Euclidean distance between them is at most 1 and 2) k-nearest neighbor graphs in which each node establishes 
(undirected) edges to the k points nearest to it. We will refer to the former model as UDG(2, A) and the latter 
as NN(2,fc) (following the notation introduced in [8]). The 2 in both these terms denotes the dimension of the 

(21 

space. Both these models display a critical phenomenon. For UDG(2, A) it is known that there is a value Ac 
such that if A > A^ 2 ' the graph contains an infinite component [TCj. For NN(2, k) the density is not relevant, 
instead Haggstrom and Meester show that k is the critical parameter (see [5]) i.e. there is a k c (2) such that for 
all k > k c (2), NN(2, k) has an infinite component. We show that as long as the critical parameters take at least 
certain values (which are higher than the critical values) it is possible in both cases to construct a subgraph of 
the infinite component with the following properties: 

PI. (Sparsity) The subgraph has a maximum degree of 4. 

P2. (Constant stretch) The distance between any two points in the subgraph is at most a constant factor greater 
than the Euclidean distance between the points. 

P3. (Coverage) The subgraph is infinite and the probability of a square region of K 2 not containing any points 
of the subgraph decays exponentially with the size of the region. 

P4. (Local computability) Each node can determine if it is part of the subgraph by using its location information 
and by communicating with its immediate neighbours. 

The subgraph constructed for UDG(2,A) is called UDG-SENS(2, A) and that constructed for NN(2, k) is called 
NN-SENS(2, k). We will show that there is a value A s (and fc s ) such that for all A > X s (resp. k > k s ), 
UDG-SENS(2, A) (resp. NN-SENS(2, k)) has properties (P1)-(P4). These properties align well with the prop- 
erties of power-efficient spanners studied in the context of ad hoc wireless networks (see [HI pp 177-178].) The 
difference being that not every point of the point process is required to be part of the network as long as the 
sensing function is satisfied (which the coverage property (P3) ensures.). 

Property (P2) is of major consequence to the power consumption of the network. This follows from the rela- 
tionship between the stretch in distance between two points and the consequent increase in the power consumed 
in communicating between them. Formally, if we consider a wireless network G formed on a set of nodes V in 
R 2 , the distance stretch, S, of a subgraph H C G is defined as 

d H (u,v) 
d = max ■— -, 

u,veV dc(u, v) 

where da(u, v) is the graph distance between u and v in G i.e. distance between u and v using the edges graph 
G and djj(u,v) is the graph distance in H. Li, Wan and Wang [Til Lemma 2] showed that given a connection 
network G and a subgraph H with distance stretch 5, the power taken to communicate between any two nodes 
is at most 5@ where fi is a parameter varying between 2 and 5. Clearly a network with property (P2) achieves 
a constant power stretch since the Euclidean distance between two points is a lower bound on the distance 
between them in both UDG(2, A) and NN(2, k). Hence we claim that our constructions are power-efficient up to 
a constant factor. Additionally we have property (P3) that guarantees coverage of the region being sensed. We 
show that the probability that a region does not contain a point of UDG-SENS(2, A) (or NN-SENS(2, k)) decays 
exponentially with the size of the region when A > X s (respectively k > k s ). In both cases the decay is sharper 
if a larger value of A is chosen. This allows us to achieve a target coverage by increasing the density to a high 
enough level. 

The basic idea behind our constructions is to couple the random graph in R 2 with a discrete site percolation 
process in Z 2 . The subgraph we construct mimics the nodes of a percolated mesh. One important byproduct of 
our constructions is that they give the best known upper bounds on the critical values for both our setting. We 
improve the best known bound of 213 (due to Teng and Yao [Uj) for the critical value of k for NN(2, k) to 188. 
We also improve the bound for the critical value of A in UDG(2, A) to 1.568. 

The construction is easy to realize using location information (which can be obtained using a GPS) and local 
computation, hence satisfying property (P4). The algorithm for routing on our subgraph constructions is based 
on a simple distributed algorithm for routing on the percolated lattice given by Angel et. al. |Tj. 

In the rest of this section we introduce some notation and definitions that will be required. We also discuss 
the various strands of research relevant to our paper. In Section[2]wc describe our constructions UDG-SENS(2, A) 
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and NN(2, k) and give lower bounds on the values A s and k s above which properties (P1)-(P4) hold. The results 
regarding stretch and coverage are detailed in [5] In Section 0] we discuss the algorithmic issues involved in 
forming our subgraphs from the underlying structure and also describe how to route packets once the structures 
are made. 

1.1 Preliminaries 

We will use the notation d(x,y) to denote the Euclidean distance between two points x,y € R 2 . In general we 
will denote the graph distance (i.e. the shortest path along the edges of the graph) of two vertices u, v of a graph 
G by dc(u, v). For the random graphs NN(2, k) and UDG(2, A) we will use the notation d n (u, v) and d u (u, v) to 
denote the graph distance between them. 

Poisson point processes Our random point sets are generated by homogenous Poisson point processes of 
intensity A in R 2 . Under this model the number of points in a region is a random variable that depends only on 
its d-dimensional volume i.e. the number of points in a bounded, measurable set A is Poisson distributed with 
mean XV (A) where V(A) is the d-dimensional volume of A. Further, the random variables associated with the 
number of points in disjoint sets are independent. 

Unit disk graphs The random graph model UDG(2, A) is defined as follows: Given a set of points S generated 
by a Poisson point process in R 2 with density A, there is an edge between points x € S and y € S if d(x, y) < 1. 

fc-nearest neighbor graphs The random graph model NN(2, k) defined as follows: Given a set of points S 
generated by a Poisson point process in R 2 there is an (undirected) edge between points x € S the k points in 
S \ {x} that are closest to x. Note that the event that two points have exactly the same distance from a point 
is a measure event, but for practical purposes we can use any tie-breaking mechanism we deem fit. 

Site percolation Consider an infinite graph defined on the vertex set Z rf with edges between points x and y 
such that ||x — y\\i — 1. Site percolation is a probabilistic process on this graph. Each point of 1 d is taken to 
be open with probability p and closed with probability 1 — p. Hence we have a sample space O = ri:rez d {0' 1}j 
individual elements of which are configurations u) = (uj(x) : x E Z d ). The product of all the measures for 
individual points forms a measure for the space of possible configurations. An edge between two open vertices 
is considered open. All other edges are considered closed. A component in which open vertices are connected 
through paths of open edges is known as an open cluster. It is known that there is a value p c such that for all 
p > p c the graph obtained has an infinite open cluster. This value is known as the critical probability. When 
p > p c then each point of Z d has some non-zero probability of being part of an infinite cluster. The reader is 
referred to [7] for a full treatment of percolation and to [4] for a recent update on some new directions in this 
area. 

When the lattice undergoes percolation, the path between two connected vertices might become long and 
tortuous. We introduce some notation for this setting. The distance between two lattice points x,y £ II 1 will be 
denoted D(x, y). When the lattice has been percoalted with probability p, the distance will be denoted D p (x, y). 
Antal and Pisztora studied this setting and proved a powerful theorem which we state here as a lemma [2] 
Theorems 1.1 and 1.2]. We adopt the restatement of Angel et. al. [U Lemma 8]. 

Lemma 1.1 |H CF F° r an V P > Pc and any x,y connected through an open path in a cube M d of the infinite 
lattice. For some p, C2 > depending only on the dimension and p and for any a > p ■ D(x,y) 

pr(D p (x,y) > a)) < e~ C2a . 

1.2 Related work 

Wireless networks Topology control in wireless networks has been studied extensively (see e.g. the surveys by 
Santi |18j and Rajaraman [17j.) Two important goals of the research in this area have been ensuring connectivity 
of all nodes and energy-efficiency 
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The approach has been to take the underlying topology as a unit disk graph [19| or a proximity graph (for 
which several proposals exist c.f. surveys cited above) on a point set and then construct some kind of spanning 
subgraph of this point set with low degree, constant stretch and the property that each node can compute its 
connections using local information. Although this line of research has the same flavour as our work, it is different 
in a fundamental way - we do not require all nodes to be connected - and so we do not survey the literature 
in detail instead referring the reader to general surveys on topology control by and and a specific survey on 
spanners by Li and Wang [15] . 

Geometric random graphs and percolation The study of random graphs obtained by applying connection 
rules on stationary point processes is known as continuum percolation. Meester and Roy's monograph on the 
subject provides an excellent view of the deep theory that has been developed around this general setting [IB] . 
UDG(2, A) is studied in [IB], where the existence and non-triviality of the critical density is demonstrated. Kong 
and Zeh [TJ] show a lower bound of 0.7698 on A c . An upper bound of 3.372 was earlier shown by Hall [5]. Hall's 
paper states an upper bound 0.843 for a model of intersecting spheres which scales by a factor of 4 using scaling 
property of coninumm percolation models |16j . 

The NN(<i, k) model was introduced by Haggstrom and Meester [8]. They showed that there was a finite 
critical value, k c (d) for all d > 2 such that an infinite cluster exists in this model. They proved that the infinite 
cluster was unique and that there was a value do such that k c (d) = 2 for all d > do- Teng and Yao gave an upper 
bound of 213 for k c (d) [21]. 

/c-nearest neighbor graphs on random point sets contained inside a finite region have been extensively studied. 
The major concern, different from ours, has been to ensure that all the points within the region are connected 
within the same cluster. Ballister, Bollobas, Sarkar and Walters [3j showed that the smallest value of k that will 
ensure connectivity lies between 0.3043 logn and 0.5139 logn, improving earlier results of Xue and Kumar [22] . 
Ballister ct. al. also studied the problem of covering the region with the discs containing the k- nearest neighbours 
of the points. We refer the reader to [3] for an interesting discussion relating this setting to earlier work by Penrose 
and others. 

2 The subgraph constructions 

In this section we describe our constructions and prove some important properties. We begin by giving a general 
overview of our technique, then move on to the specifics of the two settings. For both constructions we proceed 
by viewing M 2 as a union of a countably infinite set of square tiles. Inside each tile we look for a two kinds of 
points. The first kind is what we call a representative point. Representative points lie at the centre of the tile, 
roughly speaking. We also look for relay points, which help connect representative points to neighbouring tiles. 
Both these kinds of points have precise definitions that differ for UDG-SENS(2, A) and NN-SENS(2, fc), we will 
discuss those in Section |2~T1 and |2~21 respectively. A tile in which we find both kinds of points we call a good tile, 
other tiles arc bad. We connect representative points to four relay points, one for each neighbouring tile. Several 
points within and outside good tiles may be left unconnected. See Figure [I] for a pictorial depiction. Note that 
representative points have degree 4 and relay points have degree 2. 

The subgraph drawn by connecting representative points through relay points will be the network we will use 
for sensing. In order to prove properties (P1)-(P4) we will couple the tiling with a site percolation process in Z 2 . 
We associate each tile in R 2 with a point in Z 2 . We declare a site in Z 2 open only if the tile corresponding to it 
in M 2 is good. Hence the probability that a site is open is equal to the probability that its corresponding tile is 
good. Our definitions of representative and relay points will ensure that if two neighbouring tiles are good then 
their representative points are connected through their relay points as shown in Figure [2] This corresponds to 
the edge between two open sites being open in Z 2 (see Figure [2]) 

Since paths in Z 2 corresponds to paths between points in R 2 it follows that if percolation occurs in the Z 2 
then an infinite component must exist in the geometric random graph model as well. Hence we can conclude 
that if the probability of a tile being good exceeds the critical probability for site percolation, the geometric 
random graph model also has an infinite component in it almost surely. Let us now take a more specific look at 
the constructions for UDG(2, A) and NN(2, k). 
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Figure 1 : A portion of the tiling of ' 
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Figure 2: The part of Z 2 corresponding to the tiles shown in Figure [TJ 



2.1 Unit-disk graphs 

The internal structure of a tile for the construction of UDG-SENS(2, A) is shown in Figure [3l 

We consider square tiles of side 4/3. For the sake of exposition let us assume that the tile shown in Figure [3] 
is centred at (0,0) and it's lower left corner is (-2/3, -2/3). Within each tile we consider five disjoint regions, 
the representative region Co(t), and the relay regions Ei(t), E r (t), E t {t) and Et,(t). Co(t) is a circle of radius 1/2 
centred at the origin. The regions Ei(t),i £ {l,r,t 7 b} are intuitive to understand but slightly tricky to describe 
formally. We describe one of them, E r (t). In order to do so, let us denote by t r the tile immediately to the right 
of the tile t i.e. the tile centred at (4/3,0) with bottom left corner (2/3,-2/3). It's leftmost relay region is Ei{t r ). 
Now we define E r (t) as the part of the intersection lying wholly within t of all circles of unit radius centred at 
points in Co(t) and Ei(t r ). From this set we remove all the points of Co(t). In the figure the region is depicted 
by an ellipse, but clearly it is a less regular shape. 

We call a tile t good if each of Co(t) and Ei(t), t e {I, r, b, t} contains at least one point of the point process. 
One of the points contained in Co(t) will be the representative point for this tile, denoted rep(t). Four other 
points, one from each of the regions Ei(t),t £ {l,r,b,t} will be the relay points for this good tile. If a region 
has more than one point, the tie has to be broken. This will be done in a distributed fashion. We postpone the 
discussion of this aspect to Section [4] Note that some of the relay regions overlap and hence it may be the case 
that one point fulfils two relay functions. 

According to the program described earlier we create a bijection, 0, between the tiles in R 2 and points in 
Z 2 such that neighbouring tiles in R 2 correspond to neighbouring points in Z 2 . We couple UDG(2,A) to a site 
percolation process in Z 2 by saying that a given point x in Z 2 is open only if the tile t = </> -1 (a;) is good. Now 
we can claim that the existence of an edge in Z 2 implies the existence of a path from the representative points of 
the two tiles corresponding to the two end points of the edge. We state this formally, including an observation 
about the distance stretch between the two representative points. 

Claim 2.1 // an edge exists in the percolated mesh Z 2 between two points x and y then 
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a = 4/3 
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E b 





Figure 3: A tile t and it's 5 relevant regions. Note that the region E r is the part of the intersection lying wholly 
within t of all unit discs centred at points in Co and in the region Ei (i.e. the left relay region) of the neighbour 
tile t r . 



1. There is a path between the representative points rep{4> (x)) and rep(<f> 1 (y)) of the tiles corresponding to 
x and y in UDG(2, A) and 



2. there is a constant c u < 3 such that 
d k (repicj)- 1 (x) ) , rep{<p" 1 (y ) ) ) 



< c u ■ d{rep{<p- l {x)), repi^iy)))- (1) 



Figure 4: A path between the representative points of neighboring good tiles. 



Proof. Clearly if two neighbouring tiles t and t' are good, by the goodness condition there will be an edge 
from the representative point of one of them to a relay point in the direction of its neighbor. This relay point 
will subsequently connect to the relay point of that neighbor closest to it, which will in turn be connected to the 
representative point of the neighbour (see Figure QJ. Clearly each of the three edges on the path from rep(£) to 
rep(i') is at most 1 unit in length so c u < 3. □ 

The largest connected component formed by the representative points and relay points is UDG-SENS(2, A). 

From Claim [23] it is easy to deduce that if an infinite component exists in the site percolation setting, then 
an infinite component exists in UDG(2, A). Hence we need to determine for what values of A the site percolation 
process is supercritical. The critical probability for site percolation lies between 0.592 and 0.593 (see e.g. j!3j). 
Numerical calculations showed that the smallest value of A for which the probability of a tile being good exceeds 
0.593 is A s = 1.568. Hence for A larger than this value UDG-SENS(2, A) is infinite. 

Since this improves the best known upper bound of 3.372 [12] . we state it here as a theorem. 
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Theorem 2.2 For UDG{2, A) 

A^ 2) < 1.568. 

In Section[3]we will show that UDG-SENS(2, A) for A > A s has constant stretch. For now we move on to the 
construction of NN-SENS(2, k). 

2.2 Nearest-neighbor graphs 

The internal structure of a tile for the construction of NN-SENS(2, k) is shown in Figure [5] 
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Figure 5: A tile t and it's 9 relevant regions. Note that the region E r lies wholly within all discs of the form C x 
and C z centred at points on the boundary of the discs Co and C r . 

Let us say that the tile in Figure [5] is centred at (0,0) with bottom left corner (—5a, —5a) and top right 
corner (5a, 5a). For convenience we will refer to the tiles surrounding the tile t as, couunterclockwise starting 
from the right t r , t t , ti and We consider five circles of radius a: Co centred at (0,0), C centred at (— 4a, 0), 
C r centred at (4a, 0), Ct centred at (0,4a) and C centred at (0, — 4a). There are four other region which are 
named Ei 1 E r , E t and Ef, in the figure. E r is defined as follows. Consider the largest circle centred at any point 
in Co or C r that lies wholly within the two tiles t and t r . Two such circles, C x and C z , are depicted in Figure O 
E r is the locus of the points contained in all such circles. The regions Ei, E t and E\, are defined similary by Co 
alongwith C , Ct and C respectively and the tiles ti , tt and tb respectively. 

Now, we call tile t good if 

1. the number of points inside t is at most k/2 and 

2. the nine regions Co, C r , Ct, C, C&, E r , E t ,Ei and Ef, contain at least one point each. 

One point contained in Co will be the representative point of the tile i, denoted rep(t). A point from each of 
the other 8 regions will be relay points. If these regions contain multiple points we will have to select one from 
each and discard the rest. As in the case of UDG-SENS(2, A) we postpone the discussion of how to select this 
one point each to Section [3J 

According to the program described earlier we create a bijection, 0, between the tiles in R 2 and points in 
1? such that neighbouring tiles in R 2 correspond to neighbouring points in Z 2 . We couple NN(2,/c) to a site 
percolation process is Z 2 by saying that a given point x in Z 2 is open only if the tile t = (f>~ 1 (x) is good. Now 
we can claim that the existence of an edge in Z 2 implies the existence of a path from the representative points of 
the two tiles corresponding to the two end points of the edge. We state this formally, including an observation 
about the distance stretch between the two representative points. 

Claim 2.3 // an edge exists in the percolated mesh Z 2 between two points x and y then 
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1. There is a path between the representative points rep(<p 1 (x)) and rep(4> (y)) of the tiles corresponding to 
x and y in NN(2, k) and 

2. there is a constant c^ such that 

d k (rep((/)~ 1 (x)), rep((/)~ 1 (y))) 

< c k ■ d{rep{<t>-\x)), repty-^y))). (2) 




Figure 6: A path between the representative points of two neighboring good tiles. 

Proof. The proof of the claim is depicted in Figure [6] Clearly any circle drawn from rep(i) that stays within 
t contains all of E r in it by the definition of E r . Since there are at most k/2 points in every good tile, hence 
there is an edge from rep(i) to the point guaranteed to be contained in E r , let's call it x r . We do not make any 
claims on where the edges established by x r to its neighbours lie, observing only that any point that lies in C r 
must have an edge to x r , again by the definition of E r . However, any disc centred at a point in C r that remains 
within t and t r must contain the left disc of its neighboring tile. Hence, if t and t r are both good then a path 
from rep(i) to rep(t r ) occurs. The second part of the claim is obviously true. The constant c^. can easily be 
calculated using calculus. □ 

We define NN-SENS(2, k) as the largest connected component of the graph built on representative and relay 
points. Note that unlike UDG-SENS(2, A) there are 8 relay points within each tile here and the path between 
two representative points contains 4 relay points. Also note that only the regions Ei,E r , E t or Eb can share relay 
points. The regions C;, C r , Ct and Cb do not intersect in any way. Each of the regions may contain more than 
one point of the point process. In this case one point has to be be chosen as a representative or relay as the case 
may be. This can be easily achieved by running a simple leader election algorithm [20j between the nodes in the 
region which are all connected to each other in any case. We will discuss this issue further in Section 21 

From Claim [2~31 it is easy to deduce that if an infinite component exists in the site percolation setting, then 
an infinite component exists in NN(2, k). Hence we need to determine for what settings of our parameters a and, 
more importantly, k, the site percolation process is supercritical. The critical probability for site percolation lies 
between 0.592 and 0.593 (see e.g. [13]). Numerical calculations showed that the smallest value of k for which 
the probability of a tile being good exceeds 0.593 is 188, and the value of a for which this happens is 0.893. 

Like Haggstrom and Meester's proof for the existence of a critical value [8] and Teng and Yao's proof for the 
weaker of their two upper bounds on k c {2) [21j our proof of Theorem 12.41 proceeds by constructing a coupling 
with a site percolation process on 1? . However, our construction gives a better upper bound than Teng and 
Yao's improvement of their own result (also in [21j ) to k c (2) > 213 which uses a coupling to a mixed percolation 
process. Hence we state this result as a theorem: 

Theorem 2.4 For NN(2, k), 

k c {2) < 188. 

Having described the constructions of UDG-SENS(2, A) and NN-SENS(2,fc) and having shown that there 
are values of the critical parameters for which these constructions exist, let us now proceed to show that these 
constructions indeed have constant stretch. 
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3 Stretch and coverage 



In this section we prove that our constructions have constant stretch with high probability when the critical 
parameters have high enough values. We also show that the coverage of our constructions is very good in a 
probabilistic sense. 

3.1 Constant stretch 

The argument for constant stretch of both NN-SENS(2, k), k > k s and UDG-SENS(2, A), A > A s follow similar 
lines so we present them together. For the purposes of this section we denote the tile lengths chosen for the two 
constructions as a u (= 4/3) and ak (= 0.893). In what follows, all arguments hold for both settings except where 
explicitly noted otherwise. Also in the following we implicitly assume that k > k s and A > A s . 

Let us consider any two tiles t\ and ti whose representative points rep(ii) and rep(i2) lie UDG-SENS(2, A) 
or NN-SENS(2, k). First we relate the distance in the (unpercolatcd) lattice to the cuclidcan distance between 
these two points by observing a simple fact. 

Fact 3.1 Given the constants Cu defined in Claim [2~1\ and defined in Claim HT71 then for two tiles t\, ti 

D(*(n*(*0),*(n*fr))) < ^. ^0.^)) 

c 

where c € {c«, ci.} and a G {a u ,ak} respectively. 

Fact 13.11 along with Lemma 11.11 gives us the following theorem: 

Theorem 3.2 1. For UDG-SENS(2, A) , with A > 1.568 there are constants a and c\ depending only on A 
such that 

P(d u (x 1 y)>a-D(x 1 y)) <e-^ D< - x -y\ 

2. For NN-SENS(2 7 k), with k > 188 there are constants (3 and Ci depending only on k such that 

P(d k (x,y) >(3-D(x,y)) < e —^). 

Thcorcm l3.2l is an existential result. In Section [4] we will show how to actually find the constant stretch paths 
in a distributed way with bounded overhead. For now we proceed to show that our constructions have good 
coverage. 

3.2 Coverage 

Let us consider a square region of size I x I. Let us call this B(£). We will argue that the probability that B{£) 
contains no point of UDG-SENS(2, A) (or NN-SENS(2, k)) decays exponentially with £. As before the arguments 
here also apply to both the models. We claim the following theorem: 

Theorem 3.3 1. For A > 1.568 there are constants c^,c^ depending only on A such that 

P[\B(£) n UDG-SENS{2 1 A)| = 0] < c 4 • I 2 ■ e~ C3i . 

2. For k > 188 there are constants c$,cq depending only on k such that 

P[\B(£) n NN-SENS{2, k)\ = 0] < c 5 • £ 2 ■ e"^ 1 . 

Let us denote by TB(£) the set of tiles fully or partially contained in the B{£). Let us consider the set 
(j>(TB(l)) i.e. the set of all points in 1? which are images of the tiles in TB(£) under the mapping defined earlier. 
For B{£) n UDG-SENS(2, A) to be empty, each point of <f>{TB{£)) must be outside the infinite cluster of the 
supercritical percolation process. With this insight we now refer the reader to Theorems 8.18 and 8.21 of [7] 
dealing with the radius of finite clusters in the supercritical phase. A slight modification of the proof of Theorem 
8.21 of (due to [5]) will yield the the proof of Theorem 13.31 The details are tedious and do not add anything to 
the proof described in [7] so we omit them here. Theorem 13.31 yields the following simple corollary 
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Corollary 3.4 1. There is a constant cj such that for £ > cjlogn 

P\\B{£) n UDG-SENS(2, A)| = 0] < -. 

n 

2. There is a constant cs such that for £ > cs log n 

P\\B{£) n NN-SENS(2 1 k)\ = 0] < -. 

n 

The constants cy and c% may be larger than what the network requires for its sensing function. It seems 
intuitive that adding more nodes should decrease the values of these constants, but the statement of Thcorcm l3.3l 
does not seem to provide this. However we claim this intuition is indeed satisfied for UDG-SENS(2, A) and 
NN-SENS(2, k). In order to use Theorem 13.31 to provide a guarantee of desired coverage, we argue that the 
exponential decay indeed grows sharper when the density A increases. This is true for both UDG-SENS(2, A) 
and NN-SENS(2,fc). 

For the case of UDG-SENS(2, A) it is not hard to see, since an increase in A directly leads to an increase in 
the probability of a tile being good. This, through the coupling, increases the probability of the corresponding 
site in 1? being open. In the site percolation setting the probability of a site being part of the infinite cluster, 
9(p) is known to increase montonically with the probability p of sites being open. Our claim follows because 
the increase in 9(p) is centrally involved in the (omitted) proof of Theorem 13. 31 The claim has a slightly subtler 
provenance for NN-SENS(2, k). Essentially the argument is that if we fix some value k > fc s , increasing the 
density A allows us to use tiles of smaller side length and still achieve the desired probability (> p c ) of a tile 
being good. This implies that the number of tiles within a region increases, and hence the exponential decrease 
becomes sharper. 

4 Algorithmic issues 

We now focus on the algorithmic issues involved in building UDG-SENS(2, A) and NN(2, k) and routing packets 
in them once they are built. 

4.1 Forming the networks 

After the nodes arc laid out in their positions they have to undertake four basic steps. Firstly, they have to identify 
which tile they belong to. This involves using their location information (assumed to be of the form (x, y) E R 2 ) 
and the value of the tile width (denoted a u for UDG-SENS(2, A) and au for NN-SENS(2, k)) programmed into 
the nodes. In the second step each node determines whether it belongs to one of the special regions within the 
tile as described in Section [5] In the third step all the nodes within a region communicate to elect a leader who 
is then designated as the representative point of the tile or a relay point, as applicable. In the fourth step the 
elected points of each region form connections with the leaders of their neighbouring regions. See Figure [7] for 
a formal statement of the algorithm for building UDG-SENS(2, A). The algorithm for NN-SENS(2,fc) is very 
similar. 

The function electLeader can be realized using any distributed leader election algorithm on a complete graph 
topology since all the nodes within a region can talk to each other (see e.g. The function named connect is 
simply a handshake between the two nodes mentioned. Once the calls to this function are over the set up phase 
is completed. 

Note that the algorithm of Figure [7] will not just form the largest component but will also form other small 
components. It is possible to detect how large a given component is by attempting to send packets to distant 
nodes. The nodes of a small component can then turn themselves off if they realize they are not part of 
UDG-SENS(2, A). Detecting connectivity is an area of research in itself so we do not address the issues here, 
refering the reader to some recent work in this area [10] . 

4.2 The routing algorithm 

For routing purposes the representative points of a tile act as if they are open lattice points in Z 2 . They use 
relay points to send packets to the representative points of their neighbouring good tiles (see Figure |SJ) hence 
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Algorithm construct(UDG-SENS(2, A), a) 

1. At each node v do 

(a) Determine (location,, (a:), location^ (y)). 

(b) compute id„(ai) = location^ (x)/a u 

(c) compute id„(y) = location^ (y)/a u 

(d) compute region, u 

2. For each tile t and each region r £ {Co(t), E r (t), Ei(t), E t (t), Eb(t)} do 

(a) Build S(r,t) = {v £ t | region^, = r}. 

(b) rep(i) <- electLeader(5(Co, t)) 

(c) For r £ {£ r (t), £,(*). 3t(i), 
relay(i,r) <— electLeader(S'(r, t)) 

3. For each tile t with neighbours ti,t r ,tt, tb do 

(a) For each r £ {£ r (t), ^(t), E b (t)} 
co n n ect (rep (t), relay (t, r)). 

(b) connect(rclay(i, E r (t)), relay (t, Ei{t r )). 

(c) connect(relay(i, Ei{t)), relay(t, E r (ti)). 

(d) con nect(relay(i , £7 t (t) ) , relay(t, E 6 (i 4 )) . 

(e) connect(relay(i, Eb(t)), relay(i, E t (tb))- 

Figure 7: Building UDG-SENS(2, A). 

realizing open edges in Z 2 . With this simple idea in place, we can just plug in any algorithm which performs 
routing in the percolated mesh. 

We rely on the algorithm for efficient distributed routing in the giant component of a percolated mesh given 
by Angel et. al. [T]. Their algorithm proceeds by trying to follow a shortest path from source to destination. If 
at any point the path is broken (i.e. one of the nodes is closed) they try to find the next node along the path that 
is open by performing a distributed BFS from the current location of the packet. For our purposes we assume 
that the canonical shortest path between any two nodes (xi,yi), (£2,2/2) G Z 2 is the x — y path: the path that 
proceeds to first fix the x coordinate then the y coordinate i.e. {x\,y\) — > (x2,yi) — ► (2^2, 2/2)- We describe the 
routing algorithm in Figure [9] 

The function computeNext mentioned in Figure [9] finds the next node along the x — y path. The function 
isOpen involves checking if the next tile along the path is good or not. This can be done by asking the relevant 
relay if it has a neighbour in the target tile. The function sendTo uses the relays to pass the packet to the 
representative node of the next tile which then continues the routing process. This is the simple part of the 
algorithm. If the target tile is not good a BFS is launched to find the next good tile along the x — y path. This 
is a distributed algorithm that requires nodes to be probed as the search proceeds. Finally when it finds the 
destination it has to also report the path back to the node that launched the search. Once this path is known 
the function sendToNode sends the packet along this path to the discovered node. We refer the reader to Angel 
et. al. [1] for a proof that the expected number of probes required for this algorithm is at most a constant times 
the length of the shortest path. 

5 Conclusion 

In this paper we have shown that it is possible to construct sparse power-efficient wireless ad hoc sensor networks 
with good coverage if we are willing to accept a certain level of redundancy in the system. Ideas from percolation 
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Figure 8: The path between two representative points mimics the path in Z 2 using relay points to realize edges. 



theory have been used to demonstrate that the infinite cluster of two kinds of geometric random graphs contains 
a good subgraph with the properties that we seek and that this subgraph can be built efficiently using only local 
information. It is our conjecture that the subgraphs we build should exist whenever an infinite cluster exists in 
the geometric random graphs we study. One major direction for future research involves resolving this conjecture 
one way or the other. Even if this conjecture is not true, it should be possible to bring the values of k s and A s 

(2) 

closer to the critical values Ac and k c (2). 
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